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5^ ■ Abstract 

A non perturbative computation of the evolution of singlet parton densities without gauge-fixing requires 

. a gauge invariant gluon source operator. Within the Schrodinger Functional scheme (SF), such a source 

\^ • 

' can be defined in terms of path ordered products of gauge links, connected to the time boundaries. In 

this paper we adopt this definition and perform a one loop lattice computation of the renormalization 
constants of the twist-2 operators that correspond to the second moment of singlet parton densities. 
This calculation fixes the connection between the lattice SF scheme where a non perturbative evaluation 
of the absolute normalization of singlet parton densities can be made at low energy and the MS scheme 
where one can extract the experimental values. 
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1. Introduction 



Non perturbative calculations of parton densities are the the only possibility to fix their absolute 
normalization from first principles. Lattice simulations are suitable to this purpose and various 
estimates have been presented in the literature for the first moments of valence quark densities 
[1,2,3,4]. In general, experimental values of moments of the structure functions are obtained by 
comparing production rates with theoretical cross sections in common continuum scheme at high 
energy, where perturbation theory becomes reliable. For the non singlet case, the matching of 
low energy estimates of hadron matrix elements of Wilson operators in lattice schemes with their 
experimental values has been realized within the Schrodinger Functional scheme (SF), integrated 
by a finite size recursive method to match the large gap of the energy scales involved. A crucial 
element of the calculation was the non perturbative evolution of the Wilson operator taken in 
a matrix element with a proper SF quark source that defined the SF scheme [5] . An analogous 
calculation for the gluon density and in general for the singlet parton densities involving the 
mixing between gluon and sea quark densities has not yet been attempted, given the difficulties 
in accounting the sea quark pair creation through unquenched simulation algorithms. In this 
paper we propose a defintion in the SF scheme of a gluon source that can be used to evaluate 
the non perturbative running of the mixing renormalization matrix characterising the singlet 
evolution. Using the Wilson action and the Feynman gauge, we perform a one loop calculation 
that fixes the relation between the SF singlet scheme and the MS scheme for the second moment 
of singlet densities. Such a calculation is preliminary to the non perturbative evaluation of the 
singlet densities hadron matrix elements and to their comparison with experimental data. The 
method can be extended to higher moments where the experimental information is scarce [6] 
and even a modeste precision can help fixing the gluon density at momentum fraction greater 
than 0.5 [7]. 

The paper is organized as follows: section 2 contains basic definitions of the SF scheme, singlet 
operators and SF quark source. In section 3 a gauge invariant gluon source is introduced and 
its perturbative O(5o) expansion is performed. In section 4 we define correlation functions 
involving singlet operators, perform the one loop calculation and in section 5 we extract the one 

1 



loop renormalization constants of the singlet operators. 



2. Singlet Structure Functions 

2. 1 Schrddinger Functional 

This section is only meant to recall some basic facts that have been discussed exhaustively in 
the literature [8]. The theory is set up on a hyper-cubic euchdean lattice with spacing a and 
size T X (throughout the paper we put T = L, writing T whenever it has to be recalled the 
time character of the variable). The Schrodinger functional represents the amplitude for the 

time evolution that takes into account quantum fluctuations of a classical field configuration 
between two predetermined classical states. It takes the form of a standard functional integral 
with fixed boundary conditions. Explicitly, the link variables are chosen to be periodic in space 
and to satisfy Dirichlet boundary conditions in time, 

Ak{x)U,=o = C(x), Ak{x)\,,=T = M(x) , (2.1) 

where C(x) and C"(x) are fixed boundary fields, and P projects onto the gauge invariant content 
of C"(x) [8]. Here we choose C = C = 0, leading to the boundary conditions 

Uix,k)\,,=o = l, U{x,k)\^,=T = l; A: = 1,2,3 (2.2) 

for the lattice gauge field. The only gauge transformations on the boundary time slices that 
preserve such boundary conditions are the global ones [8] . This property will be crucial for the 
definition of the gauge invariant surface sources. The quark fields are chosen to be periodic up 
to a phase in the three space directions, 

i/jix + Lk) = e'^^i/jix), ^{x + Lk) = i^{x)e-'^'' ; A; = 1,2, 3 (2.3) 

where Ok is kept as a free parameter. One can also distribute the phase to all lattice points 
by an abelian transformation on the Fermi fields, i.e. by changing the form of the usual lattice 
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derivative in 



V^i;{x) = ^[X^U{x, n)ij{x + afi) - (2.4) 



V;V(a;) = Imx) - \-'U{x - afi, i^r'i^ix - afi)] (2.5) 



where 

^^^^iae,/L. Q^^Q^ -7r<efc<7r (2.6) 

It is useful to define also the backward derivatives on the lattice 

V^(x)V^ = i[?(x + afi)U{x, m)-'A;1 - lp{x)] (2.7) 

V'(a;)^^ = ^bP{x) - i^ix - afi)U{x - afi,i^)\f,] (2.8) 
Similarly to the gauge field, Dirichlet boundary conditions are imposed on the quark fields, 

P+V^(x)U=o = p(x), P_V'(x)U=T = p'(x), (2.9) 



and 



V'(a;)P_U,=o = p(x), V(a;)P+Uo=T = p'(x), (2.10) 

The Schrodinger functional action, including Feynman rules for the fermionic part and further 
details can be found in [8,9]. 

2.2 Singlet operators 

In the continuum, moments of singlet structure functions are related, through the operator 
product expansion, to hadronic matrix elements of two kind of twist-2, gauge invariant, local 
operators of the form 



(2.11) 



where brackets [ . . . ] mean Lorentz indices symmetrization and D ^ = — D^, with 



(2.12) 



The operators in eq. (2.11) belong in the continuum to irreducible representations of the angu- 
lar momentum. On the lattice, given the lower (hypercubic) symmetry of the Euclidean lattice 
action with respect to that of the continuum (all 4-d rotations), the identification of an irre- 
ducible representation may require some particular combination of operators. This classification 

has been discussed for example in refs. [10,11]. A subset of the basis described in [11], involving 
only spatial indices is given by 



(2.13) 



Of2 = ^tr{F[i,(x)F,2](:c)} 



(2.14) 



and 



(2.15) 



Fig. 1. Graphical representation of the products of gauge field variables contributing to the lattice field tensor, 
eq. (2.15). The point x is at the center of the diagram where all loops start and end. 
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with Qfi^{x) being the sum of the plaquette loops shown in Fig. 1, maximizing the symmetry 
of the field strength on the lattice. The usual plaquette representation of the field strength 
transforms like a reducible representation of the hypercubic group, while the definition clover- 
like of eq. (2.15) tranfsorms like an irreducible one. 

The correlation functions we want to study can be organized in a 2 x 2 matrix as follows 



where Cs are the operators and S's are the quark and gluon sources, defined in the following 
sections. Beyond the tree-level there will be in general a mixing in the flavor singlet sector 

between the quark operator (2.13) and the gluon operator (2.14). However in the quenched 
approximation (Nf = 0) one has < O-^,^ §g >= 0. In this case the matrix (2.16) becomes 
triangular and the operator O'^ does not mix with O^. 

2.3 Fermionic source 

To probe operators (2.13) and (2.14) in correlation functions, one must choose suitable sources. 
In the case of the quark source we make the same choice of the non singlet calculation [13]. 
As already stressed in that paper, the SF scheme allows us to define a gauge invariant source 
that provides a spatial direction. The operator (2.13) needs two directions and we have to give 
to the quark state a momentum different from zero in one extra direction. Moreover, using a 
particular feature of the SF, we can use the constant phase term 6 defined in (2.3), called a 
finite-size momentum. The reason for this name is that, at finite volume, this phase acts like a 
momentum probed by the local operators that we want to renormalize, but unlike the standard 
lattice momentum, it escapes the quantization rule induced by the finite volume. Its value can 
be chosen smaller than the minimum value of the standard momentum Pmin = ^ ' reducing the 
associated important lattice artefacts [13]. The quark bilinear at the boundary is given by the 
expression 




(2.16) 



<OgSg> <OgSg> 




(2.17) 
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3. Gauge invariant gluon source 



In this section wc introduce a gauge invariant gluon source which can be used for the calculation 
of the singlet correlation functions [12]. To exploit the features of the SF we recall that the 
gauge group of the SF is local in the bulk and global on the boundaries. This makes the quark 
source (2.17) gauge invariant. The gauge invariant gluon source is defined by 

Sg = S = trlTiTa} (3.1) 

where the trace is over the color indices. The big-tooth state % is defined by 

^ = l^E{nKx)-nl(x)| (3.2) 

and 

1 /T \ ° 

ni(x) = TT C/o(xo,x) C/J -,x n Uo\xo,^ + ai) (3.3) 

It is natural to define a gluon source also associated with the boundary xq = T. 

S' = tT{T{T^} (3.4) 

^/ = i7E{n:(x)-n'l(x)} (3.5) 



n^(x) = — fr UQ{xQ,yL)Ui(-T,^ n U^\xo,^ + ai) (3.6) 

Eq. (3.4) is necessary in order to define the correlation function which will be used to remove 
the additional divergences of the singlet correlation functions, introduced by the source. A 
graphical representation of the source is reported in Fig. 2. It is a product of temporal links 
in the time direction from = to = r/4, connected with a spatial link, or from = T 
to xo = 3r/4. It is gauge invariant, has two spatial directions (the same ones of the operators 



(2.13) and (2.14)), and is projected at zero momentum. The gluon source at xq = T cannot be 
obtained by substituting in S xq = T/4 with xq = 3T/4. The source gives raise to linear and 
logarithmic divergences: our calculation in perturbation theory shows that they are removed by 
a proper normalization of the correlation functions where it appears. 

I — <■ — 1 xo = T/4 



Xo 



Xo = 



X X + a« 



Fig. 2. Graphical representation of Hi (x) . 

3. 1 Perturbative expansion of the gluon source 

The perturbative expansion of the gluon source is straightforward. In the SF we adopt as usual 
the time-momentum scheme defined in [9] . For the color matrices we follow the same convention 
of [9]. So we have 



U{x,^l) = exp{agoq''^{x)T-} 



(3.7) 



(3.8) 



(3.9) 



The sum over the momenta p runs in the range — vr/a < Pk < tt/o with 



27r 



(3.10) 



The Feynman rules associated with the gluon sources are found expanding in powers of go 



(3.11) 



(3.12) 



Each term of the expansion is given by the sum of various contributions that we enumerate 
alphabetically with big case latin letters (where needed). Here are the tree-level expressions: 



(3.13) 



m _ 



where repeated indices are summed. The O{go) terms are: 

T/A-a 

^ E 

p uo=0 



(3.14) 



(3.15) 



T/4-a 



2L3 



p uo,vo=0 



(3.16) 



i{l,A) 



T-a 



-:^E« E cos(|p,) g7(^,p) gSK,-P)r'''=r^ 

P uo= 3T/4 V / V 



(3.17) 



T-a 



^'^''^^ = ^ E E P^ K, P) g1(^0, -P) /"'^r^ 



P uo,vo= 3T/4 



And the 0{gQ) are given by: 
p.q 



(3.18) 



(3.19) 



T/4-a 

P.q wo>fo=o 



cos 



a 



%"K,p)g1( j,-p-q)gS(t;o,q) T-T'T^ (3.20) 
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T/4-a 



L6 



p,q uo<vo = 

T 



'-{Pi + Qi 



(3.21) 



X 9^ 4'-P-q) +cos( ^P.] j,p)g^(^;o,q) gSK,-p-q)|>r'^r''r= 



/(2,A) 



I a \ - 

p,q 



3r 



,/3T 



4 'Pj ^UT'^j'^n 4 



3r 



p - q ) T^r'T'^ 



(3.22) 



L6 



cos 



P,q MO, 1)0= 3T/4 

arjihrric 



':iPi - Qi) 



[uo,P)qi{ ^,-p-q )x 



(3.23) 



T-a 

^/(2,C) ^ __ ^ c{uo,Vo) -[COS 



L6 



P,q uo>fo= 3T/4 



:iPi + Q'i 



7o(^o,p)QoK,q) X 



(3.24) 



X 97 ( ^, -P - q) + cosf qt (^,P] Qoivo, q) 9S(^o, -p - q) [^T'^TT^ 



where 



c{uo,vo) 



1 if uo / ■vo 

i if 140= Wo 



(3.25) 



O{g°o) 



O(g'o) 



(0) 



r' 



7- 




(2,^) 



■(2,B) 



(2,C) 



Fig. 3. The contributions to the perturbative expansion of Ti up to O((7o) 



A graphical representation of the Feynman rules for the source at xq = is reported in Fig. 3 and 
the source at xq = T has an analogous graphical interpretation. From (3.13) and (3.14) one can 
see that the tree-level of the gluon sources is a gluon field at time xq = T/4 or xq = 3T/4 with 



spatial polarization and zero momentum. We have now all the ingredients that allow to perform 
a perturbative one loop calculation of the renormalization constants of the singlet operators and 
of the correlation functions involving only the sources. 



4. Renormalization 

The renormalization condition connects the bare operators on the lattice to finite operators 
renormalized at a scale ji = 1/L : 

O^ili) = Zik{iia)Ok{a) (4.1) 

As we have discussed above, in the flavor singlet sector (and in the unquenched case) there is a 
mixing between the quark (2.13) and the gluon operator (2.14) with the same quantum numbers. 
We thus write 

Of = Z,gOg + Zg,Og (4.2) 

0^ = Z,,0, + Zg,Og (4.3) 

We adopt the following renormalization conditions in the SF scheme: 

< Ofifi) > = < ^9 > I*'"' (4-4) 

< Ofifi) > l=,/L = < 0,{a) > 1*-^ = (4.5) 

< Ofifx) > = < Og{a) > 1*-^ = (4.6) 

< Of (/x) > \,=^/r. = < 0,ia) > 1*-^ (4.7) 
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From eqs. (4.2) - (4.3) and (4.4) - (4.7), the renormalization constants can be written for 
a/L < 1 in the form [14] 



Z,,{a/L) = 1 - 



-log- + 5,, + - 



+ 0(a|) 



(4.8) 



Zqg{alL) 



47r ^ 



4 , a „ a 
-log- + S,, + - 



+ o(4) 



(4.9) 



as. 



16 , a „ ^ 



+ o(«l) 



(4.10) 



Z,,{a/L) = l-^[Nf[Uog^ + BQ+N,Blg + ol^^] + 0{a%) (4.11) 

where Cp =4/3 and as = 5'^/47r. The coefficients of the logarithms represent the anomalous 
dimensions of the corresponding operators. They are responsible for the RG-evolution of the 
coefficient functions. The B's are fixed by the renormahzation conditions of eqs. (4.4) - (4.7). 

4-1 Correlation functions 

In order to define the correlation functions for the calculation of the singlet renormalization 
constants, according to eq. (2.16), we introduce the following four correlation functions 



f„{xo,P^) = -a6^e^P-(y-^)(iv;(x)7[i5 2]V'(x)C(y)72C(z)> 



f,,ixo,P^) = -a^^e^P-(^-^)(^tr{F[i,(x)F,2i(x)}C(y)72C(z)> 



(4.12) 



(4.13) 



fgq{xo) = {-tp{x)'yiiD 2]'^P{x)tl{TiT2}) 



(4.14) 



fasi^o) = {Y,HFlip{x)Fp2]ix)}tr:{T,T2}) 



(4.15) 
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In the correlation functions (4.12) and (4.13) we perform the computation with p = and 
6 = (^1, 0, 0). Moreover, we are free to choose the physical distance xq of the operator insertions 
from the lower boundary, and we fix xo = T/2 in all the correlation functions. We have also to 
define correlation functions which involve the sources both at xq = and xq = T, because the 
correlation functions (4.12) — (4.15) have to be properly normalized by removing the renormal- 
ization of the sources. The quark source has a well known logarithmic divergence [9,13]. Our 
gluon source, as it will be seen, has a leading linear divergence. The correlation function for the 
quark source is 

,12 



= E (C'(u)75C'(v)C(y)75C(z)) 



(4.16) 



u,v,y,z 

and for the gluon source: 



(4.17) 



where S and S' are defined in (3.1) and (3.4). We can calculate analytically the tree-level of 
the diagonal correlation functions fgg and fgg and of the sources correlation functions /i and 
Gi . The tree-level of /i with p = 6 = and mo = is 



A'' = AT. 



(4.18) 



where N^. is the number of colors. The tree-level of Gi is 



(0) 



' Nl-1 
1024 



(4.19) 



The tree-level of Gi is a constant and therefore the same on the lattice and on the continuum. 
The correlation functions with a non zero tree-level are fqq{xQ;0-\_) and fgg{xo). The tree-level 
of fqq is 



and in the continuum and chiral limit it takes the form 

1 



ipo) (M_(p+)e-2'^(P^)"° -M+(p+)e-2-(P+)(2r--'') 



(4.20) 



p=0 



fiq\xo;0i) 



(1+6-2^^)2 



2^{2T-xo) 



(4.21) 
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The purely gluonic correlation function has a tree-level of the form 
_ I 

fga (^o) = 2 idii (^0 + a, T/A; 0) - dn (xq - a, r/4; 0)] x 

16a^ (4.22) 

X [^22 {xo + a, r/4; 0) - ^22 {xo - a, T/A; 0)] 

where (i^^(yo) -2^0; q) is the time-momentum gluon propagator connecting yo and time slices, 
carrying a momentum q and polarization /x [9]. Eq. (4.22) is the square of the time lattice 
derivative of the spatial gluon propagator. The fact that the spatial gluon propagator with zero 
momentum is linear in time coordinates [9], implies that with ± a > r/4, the expression of 
fgg{xo) is independent of Xq and reads 

Lattice and continuum expressions are identical. 
4-2 One loop perturbative expansion 

The only missing ingredients to perform a one loop perturbative expansion of the correlation 
functions (4.12) - (4.15) and (4.16) - (4.17) are the Feynman rules coming from the expansion 
of the operators and from the gauge part of the action. In Appendix A we give the perturbative 
expansion of the field strength F^^ to first order. The second order is rather cumbersome and 
not so instructive. In Appendix B we give the Feynman rules of the gauge part of the action. 
The expansion of /i was already done in [9], while the expansion of fqq is actually the same 
performed in the non singlet calculation [13]. In the following, we will restrict the attention to 
the new correlation functions. 

i) G, 

The correlation function involving only gluon sources, Gi, defined in (4.17), can be expanded 
in terms of the perturbative expansion of S, S' and the action. The Feynman diagrams of the 
expansion of this correlation function are given in Fig. 4. There are two checks of the calculation. 
The gluon self-energy develops quadratic divergences that cancel out by summing all the self- 
energy contributions (cfr. diagrams 6. a - 6.g in Fig. 4), and the cancellation has to take place 
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separately in the gluon sector (cfr. diagrams 6. a - 
6.f - 6.g). This is verified in our computation. 
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6.e) and in the fermion one (cfr. diagrams 



The second check is specific to our definition of the gluon source. As anticipated in section 2, 
the source develops linear divergences that must cancel once we normalize fgg with \/Gi and 
they do so. Numerical results for the ratio 

^ = r[^'^ + NfR[f^^ (4.24) 

are reported in Table 1, wher it can be seen the presence of linear divergences in R^^\ 
ii) fqg 

The correlation function fqg, concerning the gluon operator with the quark source, eq. (4.13), 
can be expressed through [9] 

[C(x)C(y)]F = P-Uo{x - ab)Six,y)Uo{y - a6)P+U,=^,=„+ 

1 „ (4-25) 

2' 



-P-'jkidk +dl)a ^5xy 



which has to be expanded to the 0{g'^) order. It also requires the tree-level of the gluon oper- 
ator in the time-momentum scheme, which can be easily computed through the field strength 
expression given in Appendix A. A picture of the Feynman diagrams for fqg is reported in Fig. 5. 
For the computation of this correlation function, as in the case of fqq [13], we use a finite-size 
momentum 6 = (0.1,0,0). 






Fig. 5. Feynman diagrams of the one loop expansion of fqg. 

The expressions of the Feynman diagrams are quite involved, not so instructive and will not be 
shown. Numerical results of fqg/fqq are reported in Table 2. 

iii) fgq 

The second non diagonal correlation function, fgq is defined through eq. (4.14). Its computation 
involves only the tree-level of the gluon source. It depends on xq but not on the momentum, 
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because the gluon source is projected to zero momentum. The expansion of the quark operator 
was aheady done in [13]. Feynman diagrams are depicted in Fig. 6. Numerical results are 
parametrized as 




l.a 1.6 

Fig. 6. Feynman diagrams of the one loop expansion of /< 



f(l) 

p(o) ^^f-^gq 



NjR[J-) (4.26) 



Jgg 

because this correlation function is proportional to the number of dynamical fermions. Note that 
the value of fqg/f<^ is available for lattices with size multiple of two while R^gq"^ is available for 
lattices with size multiple of four, because correlation functions with the gluon source involve 
T/4 as an integer parameter, while for the quark source they do not. Numerical results are 
reported in Table 2. 

iv) fgg The calculation of this gluonic correlation function, eq. (4.15), is technically difficult. 
We have to use the expansion of the gluon source, eqs. (3.13) - (3.21), and of the operator up 
to the second order. We can do here the same checks that we have done for Gi. Indeed, the 
quadratic divergences cancel out by summing all the self-energy contributions (cfr. diagrams 
5. a - 5.g in Fig. 7), and the cancellation takes place separately in the gluon sector (cfr. diagrams 
5. a - 5.e) and in the fermion one (cfr. diagrams 5.f - 5.g). Linear divergences cancel out once 
we normalize fgg with \/Gi. Numerical results are presented in Table 3 by parametrizing them 
as follows 

III = 4f + (4.27) 
Jgg 

Note that Rg^g^ has a linear divergence while Rgg^^ has only a logarithmic one. 
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parameters mo and go through the renormalized ones. In a mass independent renormahzation 
scheme, we have 

9l = 9lZg{glaiJ,) (4.28) 

iriR = mqZm{gl,ayL), nig = uiq - rric (4.29) 
We choose 

mR = (4.30) 
To the order considered, the required substitution is then given by 

9o=9l + 0{g'R) (4.31) 
With niR = 0, we have 

"Zo = mW5l + 0(4) (4.32) 

For the value of mi^'' , we take the one computed in [15] 

am(^) = -0.4342856(3) (4.33) 

We are now ready to extract the renormahzation constants of the flavor singlet operators. We 
define the normalized correlation functions: 

h,g = ^, V = % (4.34) 



n v/i 



h,, = hgg = (4.35) 



The relations between the operators and the renormahzation constants (4.2) and (4.3) for the 
correlation functions are 

h''fj= ^ Za^h^^, a,f3 = q,g (4.36) 



1=1,9 
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The generic expansion of the /, h and Z functions is 



(4.37) 



/i = MO)+5o'/i(')+0(5o') (4.38) 

Z = Z(o)+ 52^^+0(5^) (4.39) 

For hqg or hgq, h^'^^ = 0. For Zqq or Zgg, Z^^^ = 1, and for Zgg and Zgq, Z^^^ = 0. By expanding 
eqs. (4.36) to order 0{gji), we have 



(4.40) 



C(^o) = gl {/^^y + ^^^^/^S^} + Oig^k) (4.41) 

/.g) (xo, e) = gl [h^S + } + 0(5^) (4.42) 

where the amphtudes (4.40), (4.41) and (4.42) on the right hand side are to be evaluated at 
mo = 0. In order to avoid technical problems without compromising the precision, we put 
mo = 1.0 X 10~^° in our programs. 

The renormalization conditions for the /I's read 

h[q^ = with xq = L/2, Oi = 0.1, /x = 1/L (4.44) 

h[f = 0, with XQ = L/2, 9x = 0.1, ^l = 1/L (4.45) 
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h[f = 0, with xo = L/2, ijl = 1/L 



(4.46) 



K?=hfJ^ ^it^ ^o = L/2, ^i=l/L (4.47) 

Since the only dependence on the lattice spacing is through the combination a/L, the continuum 
limit is equivalent to the limit L/a oo, i.e. the number of points N ^ oo. By imposing 
the renormalization conditions (4.44) - (4.47), we obtain 



J99 



Jm 



39\l' L m ^ 2 MO) ^^-^-^^ 

J 99 '^l 



5. Analysis and results 

In order to separate the logarithmic coefficients from the finite terms in the renormalization 
constants (4.47) - (4.50), and to get rid of the lattice artefacts embedded in their own definition, 
we apply a technique based on combinations of the Z's at different values oi N = L/a [9, 18]. 

From general arguments, it is known that the dependence including lattice artefacts can be 
parametrized as follows 

°° 1 L 
Z(i)(iV) = A + B log(iV) + ^ j^{Ck + Dk log(iV)}, N=- (5.1) 

fc=i 
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where all the coefficients of the expansion depend on the details of the calculation (boundary 
conditions, operator representations, etc.), except for B, which is related to the leading anoma- 
lous dimension and is therefore scheme independent. In order to check that our Z's reproduce 
the correct logarithmic coefficients, we introduce a numerical logarithmic derivative 

A(o)(Ar) = [^(i)(7v + r?)-Z«(7V-ry)] (5.2) 

with 77 = 2 for Zgq' or Zg^ and 77 = 4 for Zgg^ or Zgl^ , depending upon the presence of fermion 
or gluon sources, which are computed respectively on lattice sizes multiple of two or four. The 
quantity A^'^^N) slowly approaches B with a rate proportional to 1/N. Infact, it can be 
expanded according to 

A(°)(iV) = S + C'i5(o)(iV) + L»i/i(o)(iV) + C'2e(°)(^^) + ^2m(°)(Af) + o(^-^^ (5.3) 

with new coefficients and the auxiliary functions 
5W(iV) = l, ^W(iV) = llog(iV) 

(5.4) 

«^°^W = ]^' mW(iV) = ^log(iV) 

In order to have a safe continuum extrapolation of (5.3), terms of order 0{1/N) should be 
absent. This can be arranged with a simple procedure. The first step consists in building the 
quantity 

(,) _ A(°) (iV)ff(°) (N + rj)- A(") {N + v)gW {N) 
^ = g(0)(N + rj)-g(0){N) ^^"^^ 

which, by the same token, can be expanded as 

A(i) {N) = B + Cig^^^ (N) + 4/^^'^ (N) + D^e^^'^ (N) + O (]^) (5-6) 

where a new set of coefficients and auxiliary functions 

9^'^ (iV) = i log (1 + ^) ; /.(^) (iV), e(^) (iV) O (i,) (5.7) 
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have been introduced. Note that, at this point, the only auxiliary function of order 0(1/N) is 
g^^\N), and it can be removed with a second subtraction step, analogous to the previous one. 
By defining 

.(2).^^ ^ A(i) (iV)g(i) (iV + r?) - AW (N + 7^)gW (N) 

^ ' gW{N + n)-gW{N) ^ ' ' 

one can easily convince himself that all the terms of order 0{1/N) have been removed, so that 

A(')(N) = B + o(^-^^ (5.9) 

and the logarithmic coefficient B can be extracted with a good precision through an ordinary fit. 
Obviously, every subtraction step cancels out part of the result, and a big precision is required 
in order to avoid rounding effects. For this reason, all the computations have been done in 
double precision, as shown in Tables 1-3. Numerical results for the anomalous dimensions are 
reported on second column of Table 4, and are compared with their analytic values, obtained 
from dimensional regularization. It has to be noted that 7gg''' = and this is verified with a 
good precision by our computation. The results for jqq are omitted, because they are the same 
ones as in the non singlet calculation, and can be found in [13]. 

The extraction of the correct logarithmic coefficients shows the consistency of the calculation. 
The determination of the finite scheme dependent terms is obtained by subtracting the loga- 
rithmic divergences with the exact coefficients. This generates the subtracted renormalization 

constants 

oo ^ 

fe=l 

that are again combined in order to suppress lattice artefacts. Following eqs. (4.8) - (4.11), we 
parametrize our continuum results as 

Zqqia/L) = l + gl(j^^log^ + B^q^) (5.11) 

Z,g{a/L) = glNf (7,, log | + S^f) (5.12) 
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Zg,{a/L) = 5|(^^^log| + i?f/) (5.13) 

Z,,ia/L) = l + 9l [Nf (7(/^) log | + [bI^^^)) + ^(f log | + (5.14) 

The values of the extrapolated -Bf J (a, 13 = q,g) are reported on first column of Table 4. Again, 
Bqq^^ is omitted, because it can be found in [13]. 

The impact of lattice artefacts on the continuum approach of the finite coefficients B^^ can 
be estimated dividing Z^^^ by the continuum fit of the finite coefficients. We introduce the 
following test functions 

'^99 y^)- [Bi[]Ue) ^^-^^^ 



\a) [BSJ]i9i) 



^99[-)= ^ (5-17) 



a 



BSF 

99 



which converge to one in the continuum limit and differ from one because of the lattice arte- 
facts. Note that Z^^^ and Z^^^ coincide for the gluon-gluon correlation, because the anomalous 
dimension is zero. In each of Figs. 13 - 16 one curve refers to the unimproved case, and the 
other one shows the effect of the improvement procedure described above. After the suitable 
combinations, the continuum extrapolation is much more safe. 



The finite part of the renormalization constants can be used to match experimental results 
extracted at high energy in a specific continuum scheme like MS and lattice results calculated 
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nonpcrturbativcly at low energy and evolved to high energy within this particular lattice scheme 
(the SF scheme). The matching coefficients are given by 



^a^" = 5f J - B^f, a,P = q,9 (5.19) 
The values of B^^^ can be found in [16] and the matching coefficients come out to be 

' [Smatch](gO ^ -0.4697(1) 
[^match](/e) ^ 0.1044(1) 

(5.20) 

^match ^ -0.04162(1) 

^match ^ -0.1195(1) . 

This is the result of the calculation. 



6. Conclusions 

This computation put the premises for a lattice non perturbative calculation of the amount 
of gluon in a hadron from first principles in the SF scheme. The new definition of a gauge 
invariant gluon source might also be used for a novel definition of a^, or for further studies of 
non perturbative aspects of the gluon propagation. Preliminary numerical simulations show the 
feasibility of our definition [to appear...]. 
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Appendix A Strength tensor perturbative expansion 



In this appendix we report the perturbative expansion of the lattice strength tensor F/j^^ in 
time-momentum representation, up to order qq. Fourier transforms are required only along the 
space directions by the absence of periodic boundary conditions in time, which characterizes 
the SF scheme. This produces an asymmetry between temporal and spatial Lorentz indices and 
makes the algebra a bit more elaborate. The perturbative expansion is defined by the formula 



F«(x) = -FW(x), VfcGN 



(A.1) 



fc=0 



The 0(5q) terms are given by 



^f^^'"'^ = Ih^^""^ [cos(^)sin(ap,)g,«(xo;p) -cos(^)sin(ap,)g«(xo;p)] (a.2) 



^^'^^""^ = ^;^^e^P-|cos(^)[gl(xo +a;p) -g~,«(xo -a;p)] + 
- isiii{apk)[qo{xQ;p) + q^{xo-a;p)\ It" 



(A.3) 



and the 0(5(0) ones are 



p,q 



- cos(^) cos(^) + cos(^ + aqj) cos(^ + apfe)J ^"(aJo; v)Qi{xo\ q) + 
+ sin(^ + ^)[sin(agfc) - sin(apfc)]^j"(xo; p)gj^(xo; q) + 
+ + ^)[sin(apj) - sin(agj)]g^(a;o; p)gfc(a;o; q)} 

25 



,apj 



aqk 



(A.4) 



^ok K-^) 4 ^6 



^ e^(p+q)-x I [cos{^ + apk) + cos(^) 



2 



X [go(a;o;p)gfc(a;o +o;q) + gS(a;o -o;p)gfc(a;o -o;q)] + 
+ [cos{apk + — ) - cos(— )] X 

X [9^(2^0; p)9fe(a;o; q) + 9S(a;o -o;p)Q'fc(a;o;q)] + 

+ ism{apk)[qo{xo;p)qoixo;ci) - q^{xo - a-p)q^{xo - a;q)] + 

+ i sin(^ + ^Mixo; p)qi{xo + a; q) - q1(xo; p)g1(xo - a; q)]} r'^T' 



(A.5) 



Appendix B Feynman rules for the action 

In this appendix we report the Feynman rules in the time-momentum representation for the 
gauge part of the action, which never appeared in the hterature. The Feynman rules for the 

fcrmion one are given in [9]. The gauge part of the action is composed by a pure gauge term, a 
ghost term and a measure term. The perturbative expansion in powers of go reads 

00 

SG[q] = ^9'oSS\q] (B.l) 
k=0 

00 

Sm[q] = J29o''S^^'H<l] (B.2) 
k=l 

00 

SFp[q, c,c] = J2 9oSPp[q, c, c] (B.3) 
fc=o 

where every term in the sums of eqs. (B.l) - (B.3) has several contributions in the time- 
momentum representation, due to the fact that we must separate spatial and temporal Lorentz 
indices. The gauge action Sq and the ghost action Spp are defined in [9]. The measure term is 
defined in any standard lattice gauge theory book (see for example [17]). 
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B.l Vertices with a non zero continuum limit 

The time-momentum representation of the three-gluons vertex separates in three contributions 



i=a,b,c 

with the exphcit expressions (Feynman graphs are reported in Fig. 

T-a . 



l,m xo = a 



a 



Ij^ilk -Pk] 



X COS 



a 



ql{xQ] -1 - m)^^(xo; 1)^ • (xo; m)/ 



abc 



(B.4) 



(B.5) 



1 1 

= ;^E" E 

l,m XQ=a 



X q^ixo; -1 - m)q'j{xo; l)q'^{xo + a; m)r'' 



(B.6) 



T-a 



cos 



l.m xo=a 



Mk -Pk] 



2a 



ql{xo + a; -1 - m) + qlixo; -1 - m) 



Qo{xo;'^)qoixo;m)f 



abc 



(B.7) 





cj b j 



c] bO 




Fig. 8. Feynman diagrams of the three contributions to the three-gluons vertex. 
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The time-momentum representation of the four-gluon has many contributions and we report 
only the ones interested in computing our correlation functions. 



(2,i)r 



(B.8) 



i=a,b,c 



with the explicit expressions (Feynman graphs are reported in Fig. 9) 



X xo=a l,m,n,p 



X <! 2g«(xo;l)g-(xo;m)gj^(xo;n)g^(xo;p) x 



cos[|(Zj +mj)] cos[|(ni - cos + m^)] cos - mj)] + (^•^) 
cos[|(/j +mj)] cos[|(pi - Hi)] - 2q^{xo;\)q''j{xo;m)qf{xo;n)q'^{xo;p) x 

cos[^{ni - Pi)] cos[|k- - I tr(r"r^T=r'^) 



T-a 



X xo=a l,m,n,p 



2gf (a;o; ^)Qi{xo; m)q^{xo; n)q^{xo; p) sin(^) sin(^) + 



+ qfixo + a; l)qi{xo; m)gg(xo; n)q^{xo; p) cos [^{rii +pi)] + 
+ qi{xo;'^)qi{xo + a]m)qo{xo-,n)qo{xo-,p) cos[^{ni + Pi)] + 
- Qtixo + a; l)q^{xo + a; m)q^{xo; n)q^{xo;p) cos [|(ni + p^)] + 
+ Qoixo; n)gi (xo + a; 1)^^ (xq + a; m)^^(a;o; p) cos [^{rii + pi)] + 



- ^qt{xo-A)qo{xo;n)q^{xo + a; ni)q^{xo; p) cos[- {m + pi)] 

X tr(r'*r^r'=r'^) 
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(B.IO) 



:'(2,c)r 



T-a 



X xo=a l,m,n,p 



gi(l+m+n+p)-x ^ 



X \qt{xo;'^)qi{xo;m)q^{xo;n)qf{xo;p) x 



-10 + ^( cos(apfc) + cos(a/fc) j - 12^cos[a(/fc + mk)] 



- 2(7j''(2;o + a; l)'7i'(a;o + a; m)(?,'=(xo + a; n)^f (xq + a; p) + 
+ Sg'f (xo; 1)^ (a^o; m)Q'^^(xo; n)gf (xo + a; p) + 

+ MiixQ] \)q\{xQ + a; m)^^ (xq + a; n)gf (xq + a; p) + 

- 12g,"(xo;l)Q-(xo;m)g^(xo + a;n)gf(xo + a;p)| x 
X i^{T<'T^T''T'^) 



+ 



1 

i 


o p 


m 


° j 




KD n 


■j 





1 

i 


o p 


m 







o n 









sr\] s^r^] 

Fig. 9. Feynman diagrams of the three contributions to the four-gluon vertex. 



Also the ghost term separates in two contributions (which are depicted in Fig. 10) 
Spp [g, c, c] = ^ S^p'p^ [g, c, c] 

and these are given by 

T-a ^ ^ _^ ^ 

SFp\q-,c,c]= +mfc)cos(-mfc)5"(xo;-l-m)g^(xo;l) x 

xa=a l,m 

X ?(xo;m)/"'"= 
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T-a 



L6 

a;o = i l,m 



go(xo;l)c'=(xo;m) + 



+ qQ{xo;'\.)c''{xo + a; m) - go(xo - a; l)c''(xo - a; m) + 



(xo - a;l)c''(xo;m) 



cabc 



(B.14) 



6 



1 



1 m >° m 1 m 

Spp' [q, c, c] s'-pp' [q, c, c] 

Fig. 10. Feynman diagrams of the two contributions to the two-ghost - one-gluon vertex. 



a b 



B.2 Vertices not existing in the continuum 
The measure vertex is given by 
N, 



12a2 



EzsEf E (g~fc(xo;l)gfcK;-l)) +'Zo(xo;lK(xo;-l) 

xo=a 1 U = l,2,3 ^ ^ 



(B.15) 



Fig. 11. Feynman diagram of the measure vertex. 

The ghost action contributes also with a vertex which has a zero continuum hmit. It is made 
by two contributions (see Fig. 12): 



Spp [q, c, c] = ^ S^pp^ [q, c, c] 



(B.16) 



i=a^b 



where 



T-a 



'S'jr^P^[i?,c,c] = —a yg ^""(ico; -1 - m - n)(/fc + + nk)hkX 



12 " ^ L9 

a;o = a l,m,n 



(B.17) 



abe fcde 



X gl(xo;l)gl(xo;m)c''(xo;n)r^V 
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T-a 



S\pp^\(l, C, C 



12" ^ L9 

a;n=a l,m,n 



m — n 



oi^o; ^)qo{xo; m)c'^(xo + a; n) + 

(B.18) 



%(a;o; l)^S(2;o; m)c (xo; n) - go(xo - a; l)^o(2;o - a; m)c (xq; n) + 



+ 9o(^o - a; i}qoixo - a; m)c'^(xo - a; n) 



afae £cde 



b k 




cO 








\ \ d 1 

Oi— 1 — m — n fli— 1 — m — n 

[q, c, c] S'j?;^' [g, c, c] 

Fig. 12. Feynman diagrams of two contributions to the two-ghost - two~gluon vertex. 
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Table 1. Ratio G^^VGi"^ 
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Table 2. Ratios /^gV/i?^ and ff^jf, 
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/f(0) 

Jgg 1 Jgg 
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Table 3. Ratio /igV/< 
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Table 4. Constant and logarithmic coefHcients as obtained numerically 
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Fig. 13. Continuum approach of [^gg']'^''- Empty dots represent the original data, while the filled ones 
result from the cleaning procedure. 
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Fig. 14. Continuum approach of [■^gg']''''^'^'' • Empty dots represent the original data, while the filled ones 
result from the cleaning procedure. 
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Fig. 15. Continuum approach of ^^^g' • Empty dots represent the original data, while the filled ones 
result from the cleaning procedure. 
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Fig. 16. Continuum approach of Zg]j' . Empty dots represent the original data, while the filled ones 
result from the cleaning procedure. 
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